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ABSTRACT 

DNA methyltransferase 1 (Dnmtl) reestablishes 
methylation of hemimethylated CpG sites generated 
during DNA replication in mammalian cells. Two 
subdomains, the proliferating cell nuclear antigen 
(PCNA)-binding domain (PBD) and the targeting 
sequence (TS) domain, target Dnmtl to the replica- 
tion sites in S phase. We aimed to dissect the details 
of the cell cycle-dependent coordinated activity of 
both domains. To that end, we combined super- 
resolution 3D-structured illumination microscopy 
and fluorescence recovery after photobleaching 
(FRAP) experiments of GFP-Dnmt1 wild type and 
mutant constructs in somatic mouse cells. To inter- 
pret the differences in FRAP kinetics, we refined 
existing data analysis and modeling approaches to 
(i) account for the heterogeneous and variable distri- 
bution of Dnmtl -binding sites in different cell cycle 
stages; (ii) allow diffusion-coupled dynamics; (iii) ac- 
commodate multiple binding classes. We find that 
transient PBD-dependent interaction directly at rep- 
lication sites is the predominant specific interaction 
in early S phase (residence time T res <10s). In late S 
phase, this binding class is taken over by a substan- 
tially stronger (T res ~22s) TS domain-dependent 
interaction at PCNA-enriched replication sites and 



at nearby pericentromeric heterochromatin sub- 
regions. We propose a two-loading-platform-model 
of additional PCNA-independent loading at post- 
replicative, heterochromatic Dnmtl target sites to 
ensure faithful maintenance of densely methylated 
genomic regions. 

INTRODUCTION 

DNA methylation is an essential epigenetic mechanism in 
mammals involved in gene regulation, genomic imprint- 
ing, X inactivation and carcinogenesis (1-3). Once estab- 
lished de novo during cell differentiation, the genomic 
methylation pattern is maintained by the DNA 
methyltransferase 1 (Dnmtl), a 183 kDa-sized enzyme 
that transfers methyl groups to hemimethylated substrate 
CpG sites generated during DNA synthesis in S phase 
(2,4,5). Hence, it seems obvious that the regulation of 
Dnmtl is tightly coupled to DNA replication. In fact, 
using fluorescence recovery after photobleaching 
(FRAP), we have previously reported that Dnmtl associ- 
ates with replication foci (RF) by a highly transient inter- 
action with the replication clamp proliferating cell nuclear 
antigen (PCNA) in early S phase via the PCNA-binding 
domain (PBD) of Dnmtl, enhancing the efficiency of 
covalent complex formation at its substrate sites (6,7). 
PCNA forms a homotrimeric ring around the DNA at 
replication forks and operates as a quasi-immobile 
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loading platform for various replication-associated factors 
(8,9). Furthermore, the targeting sequence (TS) domain 
has been found to mediate association of Dnmtl to con- 
stitutive heterochromatin from late S phase into G2 
(10,11), the latter finding challenging the strict coupling 
to the replication process. In addition, the TS domain was 
implicated in the interaction with Uhrfl (12). Uhrfl is an 
essential cofactor in the DNA methylation process and it 
has been proposed that Uhrfl targets Dnmtl to 
hemimethylated sites (13-17). Although the role and regu- 
lation of Dnmtl has been a popular field of research for 
many years, the details and functional implications of the 
cell cycle-dependent coordinated binding activity of 
the PBD and the TS domain still remain elusive. 

FRAP techniques offer an effective tool to study in vivo 
the mobility of cellular proteins and to gain a better 
understanding of molecular interactions that drive or 
limit the mobility of fluorescent fusions expressed in cells 
(18-21). By bleaching a subpopulation of fluorescent 
proteins and by analyzing the redistribution of fluores- 
cence over time, one can obtain measures of the half- 
time of recovery and the size of mobile fractions. To fur- 
thermore extract kinetic parameters from fluorescence 
recovery curves, one can describe the underlying 
dynamics of the proteins by a set of differential equations 
and apply a fitting procedure. Such kinetic modeling 
approaches can be useful to detect and quantify distinct 
dynamic populations [mobility classes (MCs)] and have 
been successfully used to quantitatively characterize diffu- 
sion and to some extent interactions inside living cells 
(22-24). Within the nucleus, the mobility of protein 
factors can be limited by binding to rather immobile struc- 
tures, most prominently to chromatin, or to stationary 
enzyme complexes, such as the replication machinery, 
transcription domains or splicing speckles (25-27). The 
majority of these interactions are surprisingly transient 
to accommodate dynamic exchange, which is pivotal to 
provide cellular plasticity and efficient responses to 
external signals (28-30). 

We aimed to extract quantitative measures of the 
binding properties of Dnmtl in vivo. Unlike for many 
other nuclear factors, Dnmtl -binding sites are 
non-homogeneously distributed with association sites con- 
stantly changing their location throughout the cell cycle. 
To characterize the contribution of the PBD and TS 
domain on these changing interactions, we analyzed the 
mobility of wild type and mutants of GFP-tagged Dnmtl. 
In particular, we interpreted the differences between the 
mutants by modeling FRAP experiments where half of 
the nucleus is bleached. As the spatial distribution of the 
different binding sites is unknown a priori, we modified 
classical reaction-diffusion models in a way that diffusion 
is simplified into a two-compartment exchange model and 
binding events are averaged over the entire half-nucleus. 
This model also takes into account potential multiple 
binding partners of proteins with different binding 
affinities. 

Our results provide evidence that the collective integrity 
of the PBD and TS domain is necessary and sufficient for 
the entire S phase-dependent targeting of Dnmtl to its 
localization sites. Transient PBD-mediated interaction at 



RF is the predominant specific interaction in early S 
phase, while in late S phase, this binding class is relegated 
by an ~2-fold stronger TS domain-dependent binding. 
Supported by super-resolution imaging with 
3D-structured illumination microscopy (3D-SIM) (31,32) 
we show that TS binding is not restricted to replication 
sites but also occurs PCNA/PBD-independently at 
postreplicative constitutive heterochromatin. We propose 
a two-loading-platform-model in which the increasing 
density of hemimethylated CpG sites in conjunction with 
increased level of heterochromatin marks at 
postreplicative heterochromatin in late S phase provides 
high-affinity binding sites for TS-mediated binding of 
Dnmtl. PCNA-independent loading downstream of repli- 
cation thus provides a mechanism to ensure maintenance 
of densely methylated heterochromatic DNA sequences. 

MATERIALS AND METHODS 

Expression constructs and cell culture 

The expression constructs, lxGFP, 2xGFP, 4xGFP, 
GFP-Dnmtl wt , GFP-Dnmtl Q162E and GFP-Dnmtl ATS 
have been described previously (6,33,34). GFP- 
Dnmtl Q162E/ATS was derived from GFP-Dnmtl ATS by 
overlap extension PCR. Mouse C2C12 myoblast cells 
were cultured in DMEM supplemented with 20% fetal 
bovine serum and 50 ug/ml gentamycine. For live cell ex- 
periments, cells were seeded in Lab-Tek chamber slides 
(Nunc) oru-slides (Ibidi), using either pools of stably 
expressing cells or transiently transfected cells. 

Creation of stably expressing cells has been described 
before (6). For transient transfections, cells were grown up 
to 30-40% confluence and transfected with TransFectin 
transfection reagent (Bio-Rad) or FuGENE HD (Roche) 
according to the manufacturer's instructions. Cells were 
then incubated overnight (TransFectin) or ~40h 
(FuGENE HD) before performing FRAP experiments (6). 

Only moderately expressing cells with unsuspicious 
morphology were chosen for further analysis. The 
overall Dnmtl level of endogenous and ectopically ex- 
pressed protein was determined for all analyzed constructs 
and cell lines by immunofluorescence labeling using a 
novel Dnmtl -specific rat monoclonal antibody 5A10. 
Quantitative analysis of labeling intensities revealed on 
average 2-fold increased protein levels compared with 
non-transfected control cells confirming no major 
overexpression of the GFP-Dnmtl fusion constructs (see 
Supplementary Methods and Supplementary Figure SI 
for details on the antibody characterization and the 
immunofluorescence assay). 

EdU pulse labeling, immunofluorescence staining and 
structured illumination microscopy 

Cells, stably expressing GFP-Dnmtl fusions were seeded 
on No. 1.5H precision coverslips (Marienfeld Superior), 
formaldehyde fixed and permeabilized with ice-cold 
methanol. For labeling of postreplicative DNA, 5 uM 
5-ethnyl-2 , -deoxyuridine (EdU) was added to the growth 
medium 60min before fixation. Endogenous PCNA 
was fluorescently labeled either with a rat monoclonal 
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antibody 16D10 (35) or a mouse monoclonal antibody 
PC 10 (Abeam) and secondary antibodies conjugated to 
Alexa Fluor 594 (Invitrogen) or CF405S (Biotium). GFP 
was postlabeled with ATT0488 conjugated GFP-Booster 
(ChromoTek). EdU was detected by Cu (I) catalyzed 
cycloaddition ('click-chemistry') of 20 uM Alexa Fluor 
594 Azide (Invitrogen) diluted in 0.1 M Tris/HCl (pH 
8.6) containing 4mM CuS0 4 and 50 mM Na-ascorbate. 
Cells were counterstained with 1 (ig/ml 4\6-diamidino-2- 
phenylindole and embedded in Vectashield (Vector 
Laboratories). 

3D-SIM was performed on a Delta Vision OMX V3 
(Applied Precision) system equipped with a 100x/1.40 
NA PlanApo oil immersion objective (Olympus), 
Cascade 11:512 EMCCD cameras (Photometries) and 
405, 488 and 593 nm diode lasers. Structured illumination 
(SI) image stacks were acquired with a z-distance of 
125 nm and with 15 raw SI images per plane (5 phases, 
3 angles). The SI raw data were then computationally 
reconstructed with channel specifically measured optical 
transfer functions using the softWoRX 4.0 software 
package (Applied Precision) to obtain a super-resolution 
image stack with a lateral (x,y) resolution of ~120nm 
and an axial (z) resolution of ~300nm (31). Images 
from the different color channels were registered with 
alignment parameter obtained from calibration measure- 
ments with 0.2 um diameter TetraSpeck beads 
(Invitrogen). 

Live cell microscopy and quantitative FRAP analysis 

Live cell imaging and FRAP experiments were typically 
performed on an \J\traVIEW VoX spinning disc micro- 
scope with integrated FRAP PhotoKinesis accessory 
(PerkinElmer) assembled to an Axio Observer Dl 
inverted stand (Zeiss) and using a 63 x/ 1.4 NA 
Plan-Apochromat oil immersion objective. The micro- 
scope was equipped with a heated environmental 
chamber set to 37°C. Fluorophores were excited with 
488 nm or 561 nm solid-state diode laser lines. Confocal 
image series were typically recorded with 14-bit image 
depth, a frame size of 256 x 256 pixels, a pixel size of 
HOnm and with time intervals of 154ms. For photo- 
bleaching experiments, the bleach regions, typically with 
a length of 8-10 um, were chosen to cover the anterior 
half of the oval-shaped nucleus. Photobleaching was 
performed using two iterations with the acousto-optical 
tunable filter (AOTF) of the 488 nm and the 514nm 
laser line set to 100% transmission. Typically, 20 
prebleach and 780 postbleach frames were recorded for 
each series. In some cases, FRAP experiments were per- 
formed on a TCS SP5 AOBS confocal laser scanning 
microscope (Leica) using comparable settings as previ- 
ously described (6). 

Data correction, normalization and quantitative evalu- 
ations were performed by automated processing with 
ImageJ (http://rsb.info.nih.gov/ij/) using a set of 
self-developed macros followed by calculations in Excel. 
Details are provided in the Supplementary Methods and 
in Supplementary Figure S4. 



Mathematical model 

The mathematical models used to statistically infer the 
kinetic parameters from corrected and normalized 
FRAP datasets are based on a compartmental approach 
and biochemical kinetic principles. The model for 
diffusion-uncoupled FRAP, i.e. for molecules that 
diffuse much more rapidly than they bind or unbind, 
has previously been described (22,36). A model for 
diffusion-coupled FRAP is developed in this work and 
illustrated in Figure 4A; a similar approach has been 
taken in (37). The model considers transitions between 
the bound and the free state of a protein with association 
rate constant k on and dissociation rate constant k u jj . As 
substantiated in the Supplementary Methods, the associ- 
ation and dissociation dynamics can be expressed in terms 
of linear ordinary differential equations (ODEs) when 
replacing k on by an effective association rate constant k* on . 
The ODEs are given below. While bound proteins remain 
fixed at the respective binding sites, free proteins diffuse 
through the nucleus, thus changing their locations. 
Movements between the bleached and the unbleached 
section are modeled with a diffusion rate constant kauf. 
Its value depends on the geometry of the cell and is not 
immediately eligible for interpretation purposes. See the 
Supplementary Methods for details on the modeling of the 
movement of proteins. 

Bleached and unbleached molecules are assumed to 
behave identically, and therefore it suffices to focus on 
one type only. Hence, let F 1 ™ , P%%,, P h h °, l " ld and P%$ d 
denote the fractions of unbleached free and bound 
proteins in the bleached and unbleached sections, 
measured with respect to all unbleached proteins in the 
nucleus. These four parameters sum up to one such that 
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The recovery curve equals F = + P h £ md \ Iff - This 
term was adjusted to the data normalization procedure 
described in the Supplementary Methods and approaches 
the value one as time progresses. 

There is possibly more than one type of binding partner 
for Dnmtl, i.e. the protein may sometimes associate to a 
partner of one type and sometimes to a partner of another 
type. These partners may differ with respect to the affinity 
of Dnmtl to enter the bound state and the mean residence 
times in this state. All binding partners with identical 
or similar kinetic properties are gathered in one MC. 
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This term seems more appropriate than classes of binding 
sites (22) because different sites with identical kinetic 
properties cannot be distinguished using FRAP data. The 
number of MCs could hence be smaller than the number of 
different binding partners. Furthermore binding-unrelated 
processes like anomalous diffusion can fall into an MC. 

Suppose there are M classes of kinetically different 
binding partners for the protein of interest, labeled with 
numbers ie{l,...,M}. For all i, define and P^f'' 

as the fractions of type-z bound proteins in the bleached 
and unbleached sections, respectively, with 



pbound,i _ pboundj _|_ pbound,i ]^q[ j. pbaundj j pbound 

fraction of type-/ bound proteins with respect to all 
bound proteins. Furthermore, denote by k* oni and k 0 ff :i 
the association and dissociation rate constants corres- 
ponding to the z'th MC. Then, the recovery is described by 
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Parameter estimation 



The mathematical model contains several unknowns: The 
model parameters k* n ; , k 0 ffj and ksff, the initial values Fq, 

Pbifi' Pbifl" 1 for the components F, P f ^ e , etc. and the 

fractions //,/, // of bleached proteins, bound proteins of 
type i, etc. Due to computational effort, parameter 
redundancies and strong correlation between some para- 
meters, it is not meaningful to statistically infer all these 
unknowns simultaneously. Instead, some values were fixed 
as follows: kdiff and fu were experimentally determined 
(see Supplementary Table S3, Figure 4B and the data nor- 
malization description above). The smallest k 0 ff value was 
set to 0.005 (see the Results section), was chosen equal 
to the first value of the FRAP curve. P^j was set equal to 
fbiFf). The association rates result from the other estimates 
as k* ni = k„ !U l,( \ - pfieey pfree statistical inference of all 
remaining variables was carried out by least squares esti- 
mation. The ODEs (l)-(3) and (4)-(7) were numerically 
solved with the Euler scheme with step length 0.03, which 
corresponded to one-fifth of the observation interval. All 
software was written in R (R Development Core Team, 
2011, R Foundation for Statistical Computing, Vienna, 
Austria). 



We estimated the model parameters for each FRAP 
curve separately and compared the estimates for curves 
from the same cell cycle phase and Dnmtl construct after- 
wards. For more details about the numerics, see the 
Supplementary Methods. 

Model choice 

In our analysis, we estimate models with different 
numbers of MCs. Because the models are nested, the in- 
clusion of more MCs always leads to a better or at least 
equally good fit. However, one may ask whether the add- 
itional computational effort for multiple MCs is worth the 
improved matching of the data. At first glance, model 
choice criteria like the Akaike information criterion 
(AIC) (38) seem appropriate. In our application, 
however, the difference in the mean squared residuals 
for different models is typically small owing to parameter 
redundancies. Because of the large number of model par- 
ameters, the AIC will often favor less MCs although the 
curvature of the recovery curves is better described by 
more complex models. For that reason, we developed a 
model selection criterion that penalizes complexity less 
rigorously and is specific to our application. Due to the 
relatively small noise in the FRAP curves (Figure 3B and 
Supplementary Figure S6), we do not expect to overfit the 
data. The criterion reads as follows. 

As explained in the Results section, the up to three 
MCs are further distinguished into one or two distinctive 
mobility classes (DMCs) and up to one catalytic mobility 
class (CMC). These have to fulfill three rules: 

(1) If a DMC or CMC is present, the fraction P bound of 
bound proteins should be above a certain threshold: 
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— Abound- 



Otherwise the DMCs and CMC are discarded, and 
we assume no MCs for this FRAP curve. 
(2) Two distinct MCs should differ substantially in their 
dissociation rates. In the model with two DMCs that 
means that one should have 
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Otherwise we assume the effective number of DMCs 
to be one. 

(3) An MC only truly contributes to the model if it 
reaches a certain size: 
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Otherwise we assume the effective number of DMCs 
to be one. 

We derive appropriate values for the above thresholds 
by cluster analysis; see the Supplementary Methods for 
details and results. For each measured curve, we now 



4864 Nucleic Acids Research, 2013, Vol. 41, No. 9 



select the model that yields the best fit. This is typically the 
model with two DMCs and one CMC, but in many cases, 
the fit of the model with one DMC and one CMC is 
equally good and hence preferred. For the chosen 
model, the original number of DMCs is replaced by the 
effective number of DMCs as determined by the above 
rules (Supplementary Figure S9). Then, for each cell 
cycle phase and protein construct, the primarily chosen 
effective number of DMCs is determined. The model 
with the according number of DMCs is chosen for this 
phase and construct. Supplementary Table S3 displays 
the mean estimates for the selected model for all FRAP 
curves. These results always assume the original number 
of DMCs and do not further reduce it to an effective 
number. 

RESULTS 

Necessity and sufficiency of the PBD and TS domain for 
S phase-specific targeting of Dnmtl 

Our aim was to analyze the S phase-dependent regulation 
of the Dnmtl-binding behavior. To that end, we 
investigated four GFP fusions: wild type Dnmtl 
(GFP-Dnmtl wt ), the full-length Dnmtl carrying a point 
mutation (GFP-Dnmtl Q162E ) within the PBD (6), a 
Dnmtl mutant carrying a deletion of a highly conserved 
part of the TS domain comprising the amino acids 459- 
501 (GFP-Dnmtl ATS ) and a Dnmtl double mutant 
containing both mutations (GFP-Dnmtl Q162E/ATS ) 
(Figure 1A). To identify different cell cycle stages, we 
first co-expressed GFP-Dnmtl constructs with PCNA 
fused to monomelic red fluorescent protein 
(RFP-PCNA) in mouse C2C12 myoblast cells and 
acquired confocal mid sections of the living cells 
(Figure IB). As previously described (6), GFP-Dnmtl wt 
co-localized with RFP-PCNA at RF in early S phase. 
Co-localization with RF was also apparent in late S 
phase, when DNA of pericentromeric heterochromatin 
(pHC) is replicated. In contrast, the double mutant 
GFP-Dnmtl^ 162E/ATS was diffusely distributed within 
nuclei throughout interphase, suggesting a deficiency to 
target RF during S phase. GFP-Dnmtl ATS still 
accumulated at RF in early and late S phase, but 
showed a slightly weaker association compared with 
GFP-Dnmtl wt , indicating the activity of PBD-mediated 
targeting in all S phase stages, independent of the 
presence of the TS domain (Figure IB and Supplementary 
Figure S2). As previously reported, GFP-Dnmtl Q162E 
showed a diffuse nuclear distribution in early S phase 
but notable association to pHC replicating in late S 
phase (6). Together, this suggests that the PBD-mediated 
interaction with PCNA is necessary for the Dnmtl local- 
ization in early S phase, but evidently not for the associ- 
ation at pHC in late S phase. 

For a more detailed view on the spatial relationships of 
wild-type and mutant Dnmtl and PCNA at RF, we used 
super-resolution 3D-SIM (31,32). Owing to the -8-fold 
improved volumetric resolution of 3D-SIM (39), we 
could clearly notice subtle variations in the Dnmtl local- 
ization that escaped detection with conventional imaging 



(Figure 2 and Supplementary Figure S3). In late S phase, 
GFP-Dnmtl wt coincides to a large extent with 
immunofluorescently labeled endogenous PCNA foci in 
locally decondensed parts of otherwise homogenously 
compacted chromocenters with some Dnmtl signal ex- 
tending slightly (by a few lOOnm) beyond the PCNA 
signal. Interestingly, RF outside of chromocenters 
showed almost no enrichment of Dnmtl (Figure 2A). In 
contrast, early S phase cells showed a more balanced 
co-localization at RF (Supplementary Figure S3A). 
Co-immunostaining of non-transfected cells with the 
Dnmtl -specific monoclonal antibody 5A10 confirmed 
the same localization characteristics for the endogenous 
Dnmtl, hence excluding potential artifacts by the 
GFP-tagging or overexpression (Supplementary Figure 
S3B). As opposed to this, GFP-Dnmtl ATS precisely 
co-localized with all PCNA marked RF inside and 
outside of chromocenters (Figure 2B). GFP-Dnmtl Q162E , 
similar to the wild type, displayed an enrichment at 
chromocenter-associated RF but also in the nearby 
regions of the chromocenters that were more compacted. 
We further noted these regions to become larger toward 
the end of late S phase, indicating that TS-mediated 
binding primarily occurs at postreplicative pHC 
(Supplementary Figure S3C). Pulse replication labeling 
with 5-ethenyl-2'-deoxyuridine and co-staining with 
PCNA confirmed the association of both, wild type and 
the PBD mutant Dnmtl to postreplicative pHC 
(Figure 2D). Hence, we conclude a strict co-localization 
of the ATS mutant with PCNA at replication sites in late 
S phase, whereas both TS domain-containing constructs 
(GFP-Dnmtl Q162E stronger than GFP-Dnmtl wt ) show a 
non-strict co-localization and a tendency to bind adjacent 
postreplicative pHC. 

To gain further knowledge about the cell cycle-depend- 
ent dynamics of Dnmtl, we compared FRAP kinetics of 
mutant proteins with those of GFP-Dnmtl wt in early S 
phase, late S phase and non-S phase cells with diffuse lo- 
calization. The latter comprises mostly Gl cells but may 
also contain a smaller subset of late G2 phase cells, ac- 
cording to the different lengths of both stages. The com- 
parison was done by half-nucleus FRAP analyses to 
quantify the strength and contribution of the PBD- and 
TS domain-mediated interactions in the distinct stages 
(Figure 3A). For a thorough quantitative evaluation of 
half-nucleus FRAP data, which preceded the application 
of the mathematical model, we developed an improved 
protocol for image registration, nuclear segmentation 
and data normalization (details described in Supple- 
mentary Methods and Supplementary Figure S4). 

Initial controls revealed that the additional expression 
of RFP-PCNA influenced the kinetics of GFP-Dnmtl wt 
(Supplementary Figure S5). Hence, to avoid any biasing 
effects, we decided not to co-express RFP-PCNA, but 
instead to collectively analyze all nuclei with diffuse 
nuclear distribution of the respective GFP fusion 
protein. Control measurements of diffusely localized 
GFP-Dnmtl Q162E in RFP-PCNA co-expressing cells 
revealed no difference between the 'early S phase' and 
'Gl/late G2' group (data not shown). We quantitatively 
analyzed half-nucleus FRAP experiments of 10-20 
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Figure 1. Domain structure and subnuclear localization of GFP-Dnmtl constructs. (A) Dnmtl consists of a large N-terminal regulatory domain 
containing PBD, TS domain and a CXXC zinc finger (ZnF) domain, and a conserved C-terminal catalytic domain. The point mutation of a highly 
conserved glutamine to glutamic acid introduced within the PBD eliminates interaction with the replication machinery (GFP-Dnmtl Q162E ). A 
deletion in the central part of the TS domain (GFP-Dnmtl AT ) was introduced to abolish interaction with constitutive heterochromatin. 
Furthermore a construct containing both mutations was generated (GFP-Dnmtl Q162E/ATS ). GFP alone is used as a control for a non-binding 
protein. (B) Spinning disk confocal mid sections of GFP-Dnmtl wild type (wt) and mutant constructs in live mouse C2C12 cells co-expressing 
RFP-PCNA to mark RF. In early S phase, GFP-Dnmtl wt accumulates at RF, whereas PCNA-binding deficient GFP-Dnmtl Q162E is diffusely 
distributed throughout the nucleus. GFP-Dnmtl a is still associated with RF, but not as prominently as GFP-Dnmtl wt . In late S phase 
GFP-Dnmtl wt , GFP-Dnmtl Q162E and GFP-Dnmtl ATS accumulate at larger RF of late replicating pHC, although with slightly less strong enrichment 
observed for both mutants. GFP-Dnmtl Q162E/ATS is distributed diffusely in the nucleus throughout interphase. Scale bar: 5 um. 



datasets for each construct and categorized cell cycle 
stage(s) (Figure 3B and Supplementary Figure S6) and 
determined half-times of the recovery (? 1/2 ) and mobile 
fractions (MF) (Supplementary Table SI) 

In accordance with our previous observations (6), GFP- 
Dnmtl wt showed a moderately reduced mobility in early S 
phase (ti/2 6.3 ± 0.3 s)comparedwithGl/lateG2phase(?i/ 2 
3.3 ± 0.1s) (Figure 3B and C). In late S phase, the recovery 
was even more reduced (tin 8-3 ± 0.6 s). Recovery kinetics 
of the GFP-Dnmtl Q162E/A ^ s double mutant comprising all 
interphase stages revealed the same fast recovery kinetics 
(h /2 3.6 ± 0.3 s) as observed for GFP-Dnmtl wt in Gl/late 
G2 phase, suggesting the complete loss of any S phase- 
specific interaction. The result implies that the collective 
integrity of PBD and the TS domain is necessary and suffi- 
cient for the entire S phase-dependent targeting of Dnmtl 
to its localization sites. 



Next, we analyzed both single mutants, GFP- 
Dnmtl Q162E and GFP-Dnmtl ATS , to dissect the specific 
role of both domains in early S phase and late S phase. 
In agreement with our previous analyses of GFP- 
Dnmtl Q162E , the kinetics measured for the pooled Gl/ 
late G2 and early S phase cells showing a diffuse distribu- 
tion, was almost identical to that of GFP-Dnmtl wt dif- 
fusely distributed only in Gl/late G2 phase (? 1/2 
3.1 ± 0.2 s versus 3.3 ± 0.1 s) consistent with the loss of 
PCNA interaction in early S phase. In late S phase, 
despite localizing similar to GFP-Dnmtl wt , the kinetics 
was slightly faster (t l/2 5.4 ± 0.4 s), indicating a contribu- 
tion of the PBD to the binding behavior of Dnmtl wt also 
in late S phase. However, the recovery was still slower as 
compared with the cells with diffuse localization, pointing 
toward an additional TS domain interaction. Comparing 
the wild-type construct with GFP-Dnmtl ATS , the mutant 
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Figure 2. Super-resolution imaging of heterochromatin association of GFP-Dnmtl constructs in late S phase. (A-C) 3D-SIM optical mid sections 
and z-projections from of C2C12 cells expressing GFP-Dnmtl wild type and mutant constructs immunostained with antibodies against endogenous 
PCNA. Profile plots were scaled between minimum and maximum intensity values for each nucleus. (A) GFP-Dnmtl™ 1 co-localizes largely but not 
strictly with PCNA inside ~200nm wide lacunas within otherwise densely packed DAPI-intense chromocenters of clustered pHC (inset al, 
arrowheads in profile plot 1 and inset a3). Anti-PCNA-labeled RF outside of chromocenters show only minor or no association of Dnmtl (inset 
a2, arrows in profile plot 2 and inset a3). (B) GFP-Dnmtl ATS strictly co-localizes with PCNA at RF inside and outside chromocenters (insets bl + b2 
and profile plots 3 + 4). An increased diffuse fraction is visible as small grainy evenly distributed nucleoplasmic background. (C) GFP-Dnmtl Q162E 
does not strictly co-localize with PCNA, but also associates with adjacent regions of pHC (arrowheads, inset c3 and profile plot 5). No association is 
detected in RF outside chromocenters (arrows, inset c3 and profile plot 6). (D) Additional replication labeling with a 60-min EdU pulse prior 
fixation. Association of GFP-Dnmtl wl and GFP-Dnmtl Q162E to chromocenter regions outside of PCNA foci is restricted to the bulk of EdU-labeled 
postreplicative chromatin (insets, arrows), while unlabeled, presumably not yet replicated chromocenter regions are still void of GFP-Dnmtl (insets, 
arrowheads). Scale bars: 5 urn and 1 um (insets). 
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Figure 3. Quantitative FRAP evaluation of GFP and GFP-Dnmtl constructs. (A) Representative time frames of exemplary half-nucleus FRAP 
series recorded with spinning disk confocal microscopy. Scale bar: 5 am. (B) Mean recovery curves displayed for all measured constructs and cell 
cycle stages. The inset illustrates cycle-dependent kinetics of GFP and GFP-Dnmtl wl alone (a), and in comparison with GFP-Dnmtl Q162E (b), 



GFP-Dnmtl ATS (c) and GFP-Dnmtl 



0162EMTS 



(d). GFP-Dnmtl w in cells with diffuse localization shows a decreased mobility compared with GFP. 



The GFP-Dnmtl wl mobility decreases stepwise in early S phase and in late S phase. The mobility of GFP-Dnmtl Q162E in Gl (late G2) and early S 
cells (diffuse nuclear localization, pooled) is almost identical to GFP-Dnmtl wl Gl (late G2) cells. In late S phase, a moderately increased mobility is 
observed for both, GFP-Dnmtl Q162E and GFP-Dnmtl ATS mutants compared with GFP-Dnmtl wt . Of note, despite comparable overall kinetics, both 
curves (dark green, dark orange) are clearly different in their shape. In early S phase, the t 1/2 of GFP-Dnmtl AT is reduced compared to 
GFP-Dnmtl wl . GFP-Dnmtl Q162E/ATS (all interphase stages, pooled) displays kinetics almost identical to GFP-Dnmtl wl in cells with diffuse 
localization. For clarity, error bars are omitted here, but shown in Supplementary Figure S6. (C) Half-times of recovery (ti/2) determined for 
each construct and distribution pattern. Error bars represent SEM. 
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showed also faster kinetics in late S phase (?i/ 2 6.0 ± 0.4 s), 
suggesting that the TS domain together with the PBD are 
necessary for the wild-type kinetics in late S phase. Besides 
the role of the TS domain in late S phase, this mutation 
also enhanced the mobility in early S phase 
(? 1/2 4.2 ± 0.3 s) in accordance with the fact that the 
early S phase pattern of RF association was less promin- 
ent (Figure IB). This result indicates that the binding by 
the PBD domain to RF is necessary, but not sufficient for 
the early S phase-specific localization of Dnmtl. 

We also noted that the mobile fractions (MF) of the 
wild-type construct within the observation time of 2 min 
dropped from around 100% in non-S phase to ~98% in 
early and late S phase. This directs to a small immobile 
fraction of covalently bound Dnmtl involved in the 
covalent complex formation during the enzymatic 
reaction (Supplementary Table SI). This observation is 
consistent with a rather slow speed of the enzyme 
reaction measured in vitro with hemimethylated substrate 
(40-43). A small immobile fraction (~1%) was also noted 
for the Dnmtl mutant construct. 

To test the general ability as well as differences in the 
efficiency of the investigated regulatory Dnmtl mutants to 
undergo covalent complex formation in vivo, we measured 
the time-dependent immobilization by FRAP on incuba- 
tion with the mechanism-based inhibitor 5-aza-2'- 
deoxycytidine (44) (see Supplementary Methods for 
details on the trapping assay). In agreement with the 
observed small immobile fractions, all mutants became 
immobilized albeit with variable efficiencies, with 
GFP-Dnmtl wt being already fully immobilized within 
30^15 min (corresponding to a trapping rate of ~3% 
min -1 ), followed by GFP-Dnmtl Q , GFP-Dnmtl ATS 
and GFP-Dnmtl Q162E/ATS (only -10% h _1 ) 
(Supplementary Figure S7). 

We conclude that the PBD and the TS domain are the 
only domains directly involved in S phase-specific target- 
ing of Dnmtl with respect to localization and kinetics and 
that both domains contribute to enhance the efficiency to 
initiate the catalytic reaction in vivo. To decipher the exact 
relationship of PBD- and TS domain-mediated binding in 
early and late S phase, however, a more sophisticated 
analysis is needed. 

Kinetic modeling of half-nucleus FRAP with multiple 
binding classes and diffusion-coupled dynamics 

We next sought to characterize the contribution of the 
PBD- and TS domain-mediated interactions in different 
stages of the cell cycle in a more precise quantitative 
manner. Hence, we utilized mathematical modeling to 
estimate the fraction of protein bound by these domains 
and the binding strength in the different cell cycle stages. 

The choice of the model was based on several consider- 
ations. First, to take account for the heterogeneous spatial 
distribution of binding sites that strongly varies in differ- 
ent cell cycle stages. Second, to correct for diffusion- 
related effects. Third, to deal with multiple potentially 
superimposing interactions, or binding classes, respect- 
ively, including a small fraction of protein covalently 
bound during the catalytic reaction (Figure 4A). Such a 



level of complexity goes beyond the assumptions of FRAP 
analysis based on present reaction-diffusion models (45- 
48). Therefore, we decided to use a compartmental 
approach with size-calibrated diffusion correction, 
suitable for experiments with half-nucleus bleaching, 
which ensure representative distribution of binding sites 
in all cell cycle stages. Each of the two halves was then 
considered as a well-mixed homogeneous interaction 
system where proteins can bind to different binding 
partners (22,36). 

Taking diffusion into account is especially important 
for assessing nuclear proteins as most of them undergo 
transient interactions in a diffusion-coupled behavior 
(18,46). This is also true for Dnmtl (Supplementary 
Figure S8). To approximate the diffusion of the protein, 
we introduced a size-dependent correction factor (k di fj) as 
a measure for the exchange of free molecules between the 
bleached and unbleached half to approximate the diffu- 
sion of the protein as similarly performed in (37) (Figure 
4A; see Materials and Methods). To estimate this 
exchange parameter, we performed FRAP calibration 
measurements of GFP monomers, dimers and tetramers, 
as they have known sizes and are presumably inert in cells 
(Figure 4B). The k di ff value corresponding to the size of 
GFP-Dnmtl was extrapolated from the k di ff values experi- 
mentally determined for the other three constructs in the 
same cell line using an exponential regression curve 
(Figure 4B and Supplementary Table S2). 

We noted that the normalized FRAP curves of 
GFP-Dnmtl constructs, in particular those from S phase 
cells, typically did not reach a straight plateau after 2-min 
observation time. Instead, they still followed a slight 
incline. As outlined above, this might be attributed to a 
small fraction of molecules actively involved in covalent 
complex formation during methyl group transfer. In vitro 
measurements have previously demonstrated a rather slow 
catalytic reaction of human DNMT1 on hemimethylated 
DNA in the range of 1-22 min per CpG (40^13). To 
account for this possibility, we added a class with 
variable fraction size but a small fixed dissociation rate 
(k 0 jj). This class is referred to as CMC. As opposed to 
that, classes with free k Q ff values will be referred to as 
DMC. The k 0 /f value for the CMC was chosen to be 
0.005 s -1 , which is equivalent to a mean residence time 
of 200 s (40). By fixing it, we avoided an additional free 
parameter in our fitting procedure. To further decrease the 
number of free parameters in our model, we also fixed the 
bleached fraction f bt to an experimentally determined 
value for each FRAP experiment (see Materials and 
Methods). Altogether we estimated the k a ff values of up 
to two DMCs, the fraction sizes of bound proteins (one 
CMC and two DMCs), and the remaining pool of free 
molecules (ff ree ). 

Although this modeling approach does not cover 
the whole details of our experimental system, it concen- 
trates on the characterization of the interactions while still 
integrating some essential information on the diffusion 
process and therefore provides a way to interpret the 
differences between the different forms of Dnmtl during 
the different phases of the cell cycle. In conclusion, 
our compartmental model uses an experimentally 
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Figure 4. Refined diffusion-coupled compartmental model for three MCs and determination of k^jr from GFP multimer measurements. (A) The 
nuclear compartment is divided into four main compartments: bleached and unbleached molecules in the bound or free state, respectively. The bound 
state can be subdivided into three compartments with specific properties. For Dnmtl, we choose two DMCs (DMC1 and DMC2, blue and green 
frame, respectively) and one CMC (red frame) with a fixed k 0 jj. Molecules bind and unbind with association and dissociation rates given by k oni and 
ko/f.i, respectively. In our refined modeling approach, the parameter //,/ is experimentally determined for each individual FRAP series. Migration of 
molecules is implemented in the model by introducing a new diffusion rate constant k di jj. This parameter corrects for the size-dependent exchange of 
the free molecules between the bleached and unbleached compartment. Parameter and variables entering the sets of differential equations are written 
in bold; predetermined/fixed values are indicated. (B) Quantitative FRAP evaluation of GFP multimers. The GFP mobility of the dimer and the 
tetramer decreases stepwise as compared with the monomer. The k dif f factor for the GFP constructs is estimated from the model with no DMC/ 
CMC. From these values, the corresponding k,gff factor for the size of GFP-Dnmtl is estimated using an exponential regression curve (inset 
diagram). 



determined cell type-specific correction for size-dependent 
diffusion effects and can extract up to three dissociation 
rates and the sizes of all bound pools and the remaining 
free pool. This way we took into account several 
interaction partners, diffusion-coupled dynamics and the 
irregular distribution of binding sites of the protein. 



Quantifying the properties of the PBD- and TS 
domain-mediated binding by FRAP modeling 

We applied the enhanced kinetic model to our FRAP data 
of GFP-Dnmtl wt and mutants. For each FRAP curve, the 
parameters of three differential equation models were 
determined using least squares estimation. These models 
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accounted for (i) no DMC/CMC, i.e. free proteins only, 
(ii) one DMC plus CMC or (iii) two DMCs plus CMC. 
The k dj ff values were chosen according to the respective 
protein sizes (Supplementary Table S2). For each individ- 
ual FRAP curve, the most appropriate model was 
determined based on the mean squared residuals and 
certain restrictions on the fraction sizes and magnitudes 
of the dissociation rates (see Materials and Methods). 
Then, for each construct and cell cycle phase, we identified 
the model that was preferentially chosen for the majority 
of datasets (Supplementary Figure S9). This model was 
then used to determine the final mean k a ff values, the re- 
ciprocal mean residence times T res =1/ k 0 jf as well as the 
corresponding fraction sizes for all measured constructs 
and cell cycle stages (Figure 5 and Supplementary Table 
S3). 

Model estimation for our reference FRAP measure- 
ments with GFP monomers and dimers provided a clear 
tendency of having no DMC (100% free fraction), 
whereas most FRAP curves of GFP-tetramers were best 
explained by a large free fraction (>91%) plus a smaller 
fraction of reduced mobility (<8%). Notably, the T res 
determined for this fraction was rather high and showed 
a large variation indicating a GFP-multimer-specific 
effect. Analyzing the GFP-Dnmtl constructs, we 
observed that a small fraction (/cmc 1-4%) was always 
assigned to the class of molecules potentially involved in 
the catalytic process (CMC). Estimation of the size of the 
CMC is numerically difficult owing to the small fraction 
size. We hence do not interpret those estimates here. 

For GFP-Dnmtl wt expressing cells in Gl/late G2, our 
model likewise estimated a large free fraction of ~80% 
and a population of ~19% with a relatively low mean 
residence time (T res ~8 s). In early S phase, the bound 
fraction doubled to —40% owing to binding to 
immobilized PCNA trimeric rings at replication forks (6, 
49). The mean residence time measured for this class was 
with T res ~10s slightly higher. The largest fraction of 56% 
was still assigned to the free pool. In late S phase cells, the 
bound fraction remained in a similar range with 48%. 
Importantly, concomitant with binding to pHC at 
chromocenters, the model identified two distinct DMCs: 
18% of the proteins were still bound with an intermediate 
strength (DMC1: T res ~10s), and an additional 28% with 
a substantially higher strength (DMC2: T res —22 s). 

An 18-22% fraction with consistent kinetics was con- 
stitutively present in all cell cycle stages of the investigated 
GFP-Dnmtl constructs. The nature of this constitutive 
MC remains unclear. In addition to the mutants described 
here, we performed FRAP analyses of a series of mutant 
constructs with deletions of potential interacting regions 
within the regulatory domain of Dnmtl, which included 
N-terminal truncations of various length (data not shown) 
and deletion within the ZnF domain (50). None of these 
mutants showed faster kinetics than GFP-Dnmtl wt in 
diffuse cells. We therefore attribute this constitutive class 
to an anomalous diffusive behavior (see discussion) and 
not to particular DNA/chromatin binding mediated by a 
specific domain. 

The modeling of FRAP data of the GFP-Dnmtl Q162E 
mutant revealed a modest reduction of the DMC1 fraction 



size as compared with the level of GFP-Dnmtl w in all 
measured phases (diffuse, late S) (DMC1: 10-14%; T res 
~9s). In late S phase cells, a second slower DMC with 
similar strength was still retained (DMC2: 23%; T res 
~19s). From these results, we conclude that Dnmtl 
binds PCNA at replication sites of early S phase cells 
with a mean residence time of — 10 s and with no more 
than 20-25% of the nuclear Dnmtl pool being involved 
in this reaction. The binding of the Dnmtl constructs with 
intact PBD to PCNA, as well as complete loss of the inter- 
action by introduction of the Q162E point mutation was 
confirmed biochemically by co-immunoprecipitation 
(Supplementary Figure S10). In late S phase, only a 
minor decrease in the mean residence times of the first 
and second DMC was observable for GFP-Dnmtl Q162E 
as compared with the wild type, indicating that the PBD 
does contribute, if only to a small extent, to the binding 
strength in late S phase. However, the overall bound 
fraction of molecules decreased compared with 
GFP-Dnmtl wt in late S phase (66% versus 52%) causing 
an overall faster FRAP kinetics. The double mutant 
GFP-Dnmtl Q162E/ATS did not establish any association 
pattern throughout interphase. In accordance, the ex- 
tracted kinetic properties were almost identical to those 
of GFP-Dnmtl wt in Gl/late G2 (19%; T res -8 s). 

Modeling of GFP-Dnmtl ATS in early S phase revealed a 
modest reduction of the fraction size and binding strength 
(DMC1: 30%; T res -8.5 s compared with 39%; T res -10 s 
in the wild type). This could argue for either a stabilization 
of the PCNA complex at the replication sites by the TS 
domain or may hint to the presence of a minor fraction of 
strong binding sites, which is too small to be identified as a 
distinct class. In late S phase, still only one DMC was 
identified for GFP-Dnmtl ATS (DMC1: 40%, T res ~9s), 
similar to GFP-Dnmtl wt in early S phase. This suggests 
a prevalence of TS domain-mediated binding over PBD 
mediated, provided that the conditions for TS binding are 
complied (i.e. high density of hemimethylated CpG sites in 
conjunction with heterochromatic marks; see discussion 
below). In this case, binding to PCNA does only seem 
to play a supportive role. This view is also in accordance 
with the finding that in late S phase, the Q162E mutation 
alone does not change the DMC1 or DMC2 substantially, 
but only leads to a moderate increase in free protein. 



DISCUSSION 

We addressed the complex problem of dissecting the cell 
cycle-dependent regulation of Dnmtl by super- resolution 
3D imaging, FRAP and kinetic modeling. Two main 
factors add to the complexity of the analysis. First, 
Dnmtl is a large enzyme with multiple regulatory 
subdomains, interaction partners and cell cycle-dependent 
regulation. Second, the distribution pattern of Dnmtl is 
highly variable throughout the cell cycle. Hence, we chose 
a global approach using half-nucleus FRAP to capture in 
all cases representative fractions of bound molecules and 
binding sites. Previous studies often used spot bleaching 
with a defined geometry that allowed the extraction of 
diffusion coefficients (24,51,52). However, such models 
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Figure 5. Parameters extracted from the kinetic modeling of GFP-Dnmtl constructs. (A) Mean residence times (T res , upper panel) and fractions of 
bound molecules (f) mu „d, lower panel) in three different classes (DMC1, DMC2, CMC) are displayed (not displayed free fractions add to 100% total 
amount). In all analyzed GFP-Dnmtl constructs, a fast population of molecules was identified with mean residence times between 6 s and 10 s 
(DMC1). The fraction of this fast population (DMC1) typically varies between 10% and 22%, but rises for GFP-Dnmtl wl in early S phase (40%) 
and for GFP-Dnmtl AT in early and late S phase (30 and 40%, respectively) due to the interaction with PCNA. A second, slower class (DMC2) was 
determined for both constructs with intact TS domain in late S phase with residence times varying from 19 to 22 s (DMC2) and an average size 
between 23 and 28%, respectively. The size of the CMC with a fixed T rcs of 200 s varies between 0 and 4.1%. Bars indicate SEM. (B) Parameters are 
sorted according to their targets PCNA and pHC and the constitutively present unspecified MC (constitutive class, CC). 



typically only included no more than one additional 
binding class, while we expected multiple interactions. 
Therefore, none of the previous models was immediately 
applicable to our case, which prompted us to devise a 
customized model. 



To eliminate a weak point of diffusion-uncoupled 
approaches, we further corrected for size-dependent diffu- 
sion using a calibration factor that was experimentally 
determined from measurements of GFP multimer 
proteins. Anomalous diffusion behavior has previously 
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been shown for GFP and dextran in the nucleus (53-55). 
This indicates that calculation of size-dependent diffusion 
differences according to the Stokes-Einstein relation might 
lead to wrong parameter estimates for proteins in the 
nucleus (56), especially when large size differences like 
between GFP and Dnmtl are taken into account. 

FRAP measurements are sensitive to experimental con- 
ditions and set-ups (57). We tried to correct for most con- 
ceivable external influences during image evaluation by 
using a tailored workflow of postprocessing steps (see 
Supplementary Materials and Methods). This involved 
image registration, constrained automated nuclear seg- 
mentation and three-step normalization/correction. 
Altogether, this allows us to compensate for lateral cell 
movement, nuclear import and export, bleaching due to 
image acquisition and variations in bleaching depth, 
which otherwise potentially affect the raw data and sub- 
sequent modeling results. Moreover, we modeled a slow 
CMC using a low fixed k 0 ff value, taking into account that 
a small fraction of molecules is likely to be involved in 
catalysis and thereby transiently immobilized by a 
covalent complex formation (40,43). Finally, we also 
reduced the number of free parameters by fixing the size 
of the bleached fraction to a value experimentally 
determined for each FRAP series. The number of DMCs 
was determined by model choice rules, which are oriented 
toward the numerical properties of the model. 

Although our model allows for three distinct MC, we 
possibly cannot estimate their number and properties 
beyond all doubt, mostly because two or more distinct 
interactions may fall into one MC (see Materials and 
Methods). If interaction strengths of multiple interactions 
are relatively close to each other, they may not be detected 
as separate classes but be captured as one with an inter- 
mediate mixed k a ff. The appearance of such parameter 
redundancies depends on the model and the values of 
the underlying parameters (e.g. 58). However, different 
dynamics can still be distinguished indirectly by a 
changed fraction size. An example for mixed interactions 
in one class is the similar kinetics of the non-specified con- 
stitutive class and of GFP-Dnmtl wt binding to PCNA. 
This rather small constitutive fraction may be attributed 
to one or more residual transient interactions. So far, we 
could not detect any specific subdomain of Dnmtl that 
would be responsible for a transient interaction through- 
out the cell cycle (data not shown). Thus, we tend to at- 
tribute this to an anomalous diffusive behavior within the 
nucleus that is identified as a pseudo-binding class. This 
may be caused by restrained accessibility of dense 
chromatin domains and transient trapping inside of 
small chromatin lacunas ('corralling') generating a 
'pseudo 1 binding effect and/or by unspecific transient 
binding with a broad distribution of binding affinities 
(53,54,59). Dnmtl could also be constitutively present in 
a free diffusing complex including interacting proteins like 
PCNA or Uhrfl. In fact, Dnmtl interactions have been 
described for a variety of proteins including other DNA 
methyltransferases, chromatin modifiers and transcrip- 
tional regulators (60). Interactions with high molecular 
weight complexes could potentially slow down diffusion 



of GFP-Dnmtl and thus contribute to the observed 
dynamics. 

Another limitation of our method is the precision limit 
imposed by a still large cell-to-cell variability due to (i) 
technical reasons like residual uncorrected cell motion or 
z-drift, irradiation/transfection-induced DNA damage or 
cell cycle arrest or (ii) biological reasons such as variations 
in endogenous expression and methylation levels or local 
environment. Therefore, for example, the quantification of 
the apparently small CMC gives no robust results. 
However, by estimating the kinetic parameters for each 
FRAP curve separately, we take into account this extrinsic 
noise and quantify it through standard errors. Finally, 
using experimental FRAP data, it can never be ruled out 
that the kinetics are influenced by variations in k on rates. 
However, the method does not allow assessing changes in 
the accessibility of the binding sites. These technical limi- 
tations could only be solved using large-scale simulations 
and even more complex models. 

Despite these potential shortcomings, by application of 
our method, one can still obtain a detailed picture of the 
distinct cell cycle-dependent dynamics of proteins. We 
have shown that the PBD and the TS domain are the 
only domains involved in direct S phase-dependent target- 
ing of Dnmtl and responsible for delaying its mobility. 
Furthermore, we discriminated two different MCs that 
could be matched to these two different domains of 
Dnmtl. In this study, we quantified the time they bind 
on average and found the binding via the TS domain to 
be >2-fold stronger than via the PBD, whereas the cor- 
responding fractions of bound protein were in a similar 
size range between 20 and 30%. 

In accordance with previous studies, we show that the 
more transient interaction with PCNA increases the con- 
centration of Dnmtl at replication sites to enhance the 
efficiency of maintenance DNA methylation (7). In 
addition, we have characterized the stronger binding 
properties of the TS domain. The related MC was only 
present in late S phase, when pericentric heterochromatin 
(pHC) clustered in DAPI dense chromocenters is 
replicated, suggesting a switch between PBD-mediated 
binding in early S phase to the TS domain-mediated 
binding in late S phase. The analysis of the single 
mutants, however, hinted at a somewhat more complex 
situation, as the deletion within the TS domain also 
influenced GFP-Dnmtl kinetics in early S phase and the 
mutation in the PBD influenced the GFP-Dnmtl localiza- 
tion in late S phase. This rather argues for a more subtle 
continuous change in binding balance instead of a simple 
on/off switching. Hence, association via the TS domain 
might occur also in early S phase, but at much lower abun- 
dance. Indeed a substantial minority (8/21) of early S 
phase cells could be better fitted with two respective 
distinct MCs, indicating some cell-to-cell variability, 
possibly in transition to mid S phase. On the other 
hand, PBD-mediated co-localization with PCNA is also 
observed in late S phase. However, as the fraction size 
with the respective faster off-rate is reduced to non-S 
phase level, this Dnmtl fraction may be handed over to 
form a late S phase-specific, more stable complex, such 
that the TS domain-specific off-rate becomes 
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Figure 6. Two-loading-platform-model for the cell cycle-dependent targeting of Dnmtl to RF and pHC. Schematic representation not drawn to 
scale. Closed and open lollipops indicate methylated and non-methylated CpG sites, respectively; flags indicate heterochromatin specific marker (e.g. 
H3K9me3). Dnmtl is depicted in green. The model postulates two auxiliary factors that act as immobilizing platforms under certain conditions: 
PCNA (red) that assembles as trimeric ring at the replication fork throughout S phase, and a second unspecified factor (e.g. Uhrfl, blue) that binds 
strongly to hemimethylated postreplicative heterochromatin. Independent of replication, most Dnmtl molecules (~80% in Gl, >50% in S phase) are 
freely roaming the nucleoplasm (left column). In addition, a non-specified MC with a pseudo k 0 ffj is constitutively present throughout interphase, 
which may be attributed to either non-specific binding to chromatin or transient trapping ('corralling') of the large enzyme in the nucleoplasmic 
environment. When replicating euchromatic sequences in early S phase (upper row) an additional ~20% fraction of the Dnmtl pool transiently binds 
via the PBD to immobilized PCNA rings (red donut) with a mean residence time (l/k OJ a-j) of ~10s (1). Targeting to PCNA at RF enhances the 
efficiency of a small fraction of Dnmtl to form metastable covalent complexes (k 0 ff,i) with hemimethylated CpG substrate sites in close vicinity. This 
may occur on already assembled nucleosomes, likely involving complex formation with one or several auxiliary factors (2 a), or directly on the naked 
DNA substrate adjacent to PCNA (2 b) or to nucleosomes (2 c). In late S phase, replication through chromatin with now abundant heterochromatic 
marks in conjunction with dense CpG methylation triggers the generation of high-affinity binding sites for an auxiliary protein. These may then act 
as second loading platform (dark blue pentagons) for TS-mediated binding with mean residence time (1 / k 0 ff j2 ) of ~22s involving ~25% of the 
Dnmtl pool. Formation of this transient complex with subsequent substrate binding of a small subset of molecules occurs either directly at the 
replication fork promoted by PBD-mediated targeting, or PCNA independently at already displaced postreplicative heterochromatin chromatin that 
may have escaped loading in the first instance (3). This second PCNA-independent loading complex may be assembled well into G2 phase, until all 
hemimethylated Dnmtl target sites are fully methylated, which finally triggers disassembly of the loading complex and dissociation of Dnmtl (4). Of 
note, this conceptual model is based on the differential availability of binding sites and the free interplay of forces. While higher affinity binding sites 
are occasionally generated also in early S phase, they may be too sparse to constitute a separate MC. 



predominant. In this stage, PCNA binding would thus be 
an auxiliary factor for TS-mediated binding of a Dnmtl 
subfraction. Besides, supported by super-resolution mi- 
croscopy, we provide evidence for additional PCNA/ 
PBD-independent binding to pHC. The latter seems to 
be dependent not only on the heterochromatic context 
alone, but also on the presence of hemimethylated 
postreplicative DNA as precondition. Accordingly we 
did not observe any pHC association in early S phase. It 
is tempting to speculate that the observed kinetics reflect 
the binding of the TS domain to Uhrfl, an essential epi- 
genetic factor that has previously been shown to target 
Dnmtl to hemimethylated CpG sites and to bind 
trimethylated H3K9 (13-17). However, we cannot rule 
out other/additional modes of binding of the TS domain 



to pHC. In support of a role of Uhrfl as a docking 
platform for Dnmtl, a previous FRAP study 
demonstrated a much slower recovery of GFP-Uhrfl 
compared to GFP-Dnmtl wt in mouse embryonic stem 
cells (61). 

In light of our data we propose a conceptual 
two-loading-platform model (outlined in Figure 6). 
According to this, the kinetic balance would shift from 
predominant PCNA/PBD binding in early S phase, 
toward TS-mediated binding in later S phase stages 
when replicating densely methylated heterochromatic se- 
quences. This shift would be triggered by the strongly 
increased appearance of hemimethylated CpG sites in con- 
junction with heterochromatic marks (e.g., H3K9me3). 
These would then offer the target for the formation of a 



4874 Nucleic Acids Research, 2013, Vol. 41, No. 9 



stable complex (involving e.g. Uhrfl) that acts as a second 
Dnmtl loading platform on postreplicative chromatin 
sites. PCNA-independent loading complexes may persist 
also beyond S phase, until all hemimethylated Dnmtl 
target sites are fully methylated, which in turn triggers 
complex disassembly and gradual loss of TS-mediated 
binding in G2 phase. Such a mechanism would thus safe- 
guard faithful maintenance of dense methylation at con- 
stitutive heterochromatin important for genome stability 
(3), against the backdrop of a rather slow and inefficient 
catalytic reaction (40). 

While we favor a model of free interplay of forces in 
conjunction with a cell cycle dependent varying abundance 
of high affinity binding sites, we cannot rule out an effect 
by an induced conformational change of the Dnmtl 
protein to expose the TS domain at the onset of late S 
phase. In this context, several modifications have been 
reported like acetylation, ubiquitination, phosphorylation, 
methylation and sumoylation, resulting in a change 
in activity and/or abundance of Dnmtl (60,62-64). 
For example, it has been shown that Uhrfl ubiquitinates 
Dnmtl at the C-terminal part of the TS domain (33). 
Further studies will have to address the exact interplay of 
Uhrfl and Dnmtl as a function of variable (hemi)methy- 
lation density and the role of posttranslational modifica- 
tions of Dnmtl. 

In the present study, we have reached substantial im- 
provements on the experimental conditions and workflow 
for the quantitative and qualitative evaluation of half- 
nucleus FRAP experiments. Still, extracting definite 
answers from modeling of such FRAP data remains diffi- 
cult and to some extent limited. Of note, our analysis 
involved multiple decisions on the data normalization, 
fixation of parameters, model choice, etc. Although all 
steps have been carried out with greatest care, this deter- 
ministic approach will still fall to some extent short. New 
stochastic modeling approaches may be able to realistic- 
ally take into account random events and may hence 
better explain intrinsic variability of the FRAP curves. 
Nonetheless, this article provides a framework for the 
global assessment and quantitative measurement of 
diffusion-coupled nuclear protein dynamics with heteroge- 
neous and variable distribution of binding sites, e.g. 
during cell cycle and development. 

Our approach provided new insights into the complex 
cell cycle dependent regulation of the multi-domain 
protein Dnmtl in the epigenetic network. We arrived at 
a probabilistic two-loading-platform model that provides 
a possible explanation how PBD and TS domain act co- 
operatively to faithfully maintain genomic methylation 
patterns through cell cycle and cell divisions. Further 
studies will address the mechanistic nature of the 
complex formation involving the TS domain and the tar- 
geting of Dnmtl to hemimethylated sites. 
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